Modelling spin qubits in carbon peapods 



Line; Barbara Montanari. 2 John H. Jefferson, 3 David G. 

Pettifor, 1 Nicholas M. Harrison, 2,4 and G. Andrew D. Briggs 1 

'Department of Materials, University of Oxford, Oxford 0X1 3PH, United Kingdom 
STFC Rutherford Appleton Laboratory, Didcot, Oxfordshire 0X11 OQX, United Kingdom 
3 Sensors and Electronics Division, QinetiQ, St Andrews Road, Malvern, Worcs. WR14 3PS, United Kingdom 
Jf Department of Chemistry, Imperial College, London SW7 2AZ, United Kingdom 

(Dated: February 2, 2008) 

We have calculated electron spin interactions in chains of Sc@Cs2 endohedral fullerenes in isolation 
and inserted into a semiconducting or metallic single-walled carbon nanotube to form a peapod. 
Using hybrid density functional theory (DFT) , we find that the spin resides mainly on the fullerene 
cage, whether or not the fullerenes are in a nanotube. The spin interactions decay exponentially 
with fullerene separation, and the system can be described by a simple antiferromagnetic Heisenberg 
spin chain. A generalised Hubbard- Anderson model gives an exchange parameter J and a Coulomb 
parameter U in good agreement with the DFT values. Within the accuracy of the calculations, 
neither semiconducting nor metallic nanotubes affect the interactions between the fullerene electron 
spins. 

PACS numbers: 73.61. Wp, 73.63.Fg, 71.15.Mb, 71.23.An, 03.67.-a 



Spin qubits have potential for controlled interac- 
tions [l[ for quantum computing. Carbon is a candidate 
host for spin qubits because in 12 C materials the small 
spin-orbit coupling and absence of hyperfine coupling en- 
sures long spin coherence times. Carbon peapods Q, 
that is, single- walled carbon nanotubes (SWNT) con- 
taining fullerenes, have been proposed as hosts for spin- 
qubits 0]. The fabrication of nanoscalc electronic de- 
vices such as field effect transistors with carbon peapods 
containing various endohedral fullerenes is well estab- 
lished When spin active metallic atoms such as Sc 
are incarcerated in a carbon cage, the system develops 
hybridized orbitals resulting in an unpaired electron de- 
localized across the fullerene cage Here we report on 
detailed numerical simulations which establish the nature 
of the spin-spin interactions both between endohedral 
fullerenes and between fullerenes and nanotubes. The 
dominant interaction is of the Heisenberg form, which 
is known to be suitable for quantum computing (QC) 
in one-dimensional chains [f|. Moreover, these results 
support this system to be promising for experimentally- 
proven QC protocols which allow chains formed of iden- 
tical units to be controlled "globally" 0], i.e. without 
the need to target individual units, leading to a scalable 
molecular quantum computer. 

Four fundamental problems need to be understood 
in order to demonstrate well-defined qubits in carbon 
peapods as proposed in Ref. 0: (i) the charge arrange- 
ment within the carbon peapods; (ii) the electron spin 
distribution; (iii) the coupling between spin-qubits; (iv) 
the nature of the spin interactions between fullerenes 
and nanotube. Density functional theory (DFT) cal- 
culations have been reported on several model systems, 
such as SWNT containing C fi n, K. T Cfin, Y@Ceo, Cs2 and 
La@C 82 d M, ELI EM 111- To date little work 
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FIG. 1: (Color online) (a) Eigenspectrum of relaxed 
Sc@C82 molecule. Left (right) hand side refers to spin-up 
(-down) electrons. Solid (dashed) lines refer to occupied 
(unoccupied) energy levels, (b) HOMO of relaxed Sc@C82- 
(c) (LUMO+f) of relaxed Cs2. Dark grey (red online) and 
light grey (blue online) lobes represent positive and negative 
phases. 



has been done on understanding the spin properties of 
peapods and the above mentioned issues (ii) , (iii) and (iv) 
have not been addressed in detail thus far. We have com- 
puted the charge and spin distributions and the electronic 
structures within Sc@Cs2 chains and peapod structures. 
We find well-defined spin-1/2 qubits on the fullerenes, 
with strong evidence for a nearest-neighbor Heisenberg 
exchange interaction. We can tune this exchange inter- 
action strength by controlling the inter-fullerene spacing 
during synthesis. In order to describe the influence on the 
spin-qubits localized on the fullerenes from propagating 
electrons or holes in the nanotube, we need to go beyond 
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DFT to a model which is capable of describing the low- 
energy charge-spin excitations of the system. We con- 
jecture a generic Hubbard- Anderson model, which has 
these properties, to estimate the low-energy spin inter- 
actions; in particular the Heisenberg exchange between 
spins along the fullerene chain and the Kondo exchange 
interaction between localized spins on the fullerenes and 
spins of propagating electrons or holes in the nanotube. 

As model systems, we choose Sc@C 8 2 in (14,7) 
(semiconducting) and (11,11) (metallic) SWNTs. Our 
calculations predict an exothermic encapsulation of 
Sc@Cs2 for both tubes. The repeat units containing one 
Sc@Cg2 molecule in the (14,7) and (11,11) peapods are 
11.42 A and 12.47 A, respectively 0. In order to study 
the exchange interaction, we use double unit cells con- 
taining two Sc@C 8 2 molecules. The interwall separations 
between the (14,7) and (11,11) tubes and Sc@C 8 2 are 
3.35 A (van dcr Waals distance) and 3.55 A, respec- 
tively. The relaxed structure of Sc@Cg2 is found to be in 
agreement with Ref. [B| with a Sc-C distance of 2.26 A. 
The DFT calculations are performed with the hybrid ex- 
change density functional B3LYP [if 
mented in the CRYSTAL package [li 
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The calculations 
reported here are all-electron, i.e., with no shape approx- 
imation to the ionic potential or electron charge density. 
The geometry optimizations are performed using the al- 
gorithm proposed by Schlcgcl et al [2(|. The crystalline 
wave functions are expanded in Gaussian basis sets of 
double valence quality (6-21G* for C and 864-11G* for 
Sc). Atomic charges are estimated using Mullikcn pop- 
ulation analysis [2l|. The system is modelled as a one 
dimensionally periodic array and reciprocal space sam- 
pling is performed on a Monkhorst-Pack grid containing 
30 symmetry irreducible k-points which converges the to- 
tal energy to within 10 -4 eV per unit cell. 

FigQ](a) shows the calculated electronic eigenspectrum 
for the relaxed Sc@C 8 2 molecule. Sc has three valence 
electrons and the ground state of Sc@C 8 2 is found to be 
a spin- 1/2 system. The unpaired electron occupies the 
highest occupied molecular orbital (HOMO) of Sc@C 8 2, 
which constitutes the spin qubit. The HOMO-LUMO 
gap is 0.530 eV. The separation between the HOMO and 
(HOMO-1) is 0.789 eV. Therefore, the HOMO is well sep- 
arated from the energy levels above and below leading to 
a well-defined qubit. The Sc@C 82 HOMO is delocalized 
across the fullerene cage as depicted in FigQ] (b), which 
is in agreement with the results found in Ref. @, I22I ]. 
Furthermore, we establish that the HOMO of Sc@C 82 is 
virtually identical to the (LUMO+1) of C 8 2, as shown 
in FigfT] (c) , whereas the lower lying orbitals are hybrids 
of Sc and C 8 2. Thus Sc acts as a perfect donor to the 
Cs2 cage for the HOMO state. 

Fig. [2] shows the electronic charge rearrangement fol- 
lowing the Sc@C 8 2 encapsulation in the (14,7) nan- 
otube. The charge depletion from the nanotube is con- 
centrated around the fullerene sites. Similar qualita- 





FIG. 2: (Color online) Charge transfer in Sc@Cg 2 @(14,7) pea- 
pod. Values for dark grey (red online) and light grey (green 
online) surfaces are ±0.001 e/A 3 . Left (right) hand shows 
front (side) views. The atom colored in light grey (gold on- 
line) is Sc. 



tive results arc obtained for the (11,11) nanotube, re- 
sembling that in La@C 82 @(17,0) @]. Table 1 shows the 
charge and spin populations in Sc@C 8 2, Sc@C 8 2@(14,7) 
and Sc@C 8 2@(ll,ll) peapods. In Sc@C 82 , 1.64 elec- 
trons transfer from the Sc atom to the C 8 2 cage in- 
dicating a partially covalent Sc-cage bond. In both 
Sc@C 82 @(14,7) and Sc@C 82 @(ll,ll) peapods, electron 
transfer occurs from the nanotube and the Sc atom to 
the C 8 2 cage due to hybridization between the occupied 
states of the nanotube and fullerenes. The charge trans- 
fer from the Sc atom to the C 82 cage in Sc@C 82 @(14,7) 
and Sc@C 82 @(ll,ll) is very similar to that of Sc@C 82 . 
The electron spin distribution is also similar as 97% of 
the density resides on the C 8 2 cage and only 3% on the 
Sc atom. The shape of the spin density distribution 
in the peapods closely resembles that of the HOMO of 
Sc@C 8 2 illustrated in Fig. [1] (b). The charge transfer 
and spin distribution are insensitive to encapsulation. 

Table 1. Charge and spin populations in Sc@C 8 2, 
Sc@C 82 @(14,7) and Sc@C 82 @(ll,ll) peapods 
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In the predicted ground state configuration of 
the peapods, the spin direction alternates along the 
Sc@C 8 2 chain; the corresponding configuration with par- 
allel spins being higher in energy. We denote these 
configurations as antiferromagnetic (AF) and ferromag- 
netic (FM) states, respectively. Both states are found to 
be Mott insulators and have a total energy lower than 
the nonmagnetic state. For the Sc@C 8 2@(14,7) peapod 
where the inter-fullerene spacing is 3.42 A, the energy 
difference between the nonmagnetic and FM states is 
0.108 eV/cell. The exchange parameter J, defined as the 
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sured lattice spacing. 
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3: (Color online) Exchange interaction strength, J, 
as a function of inter-fullerene separation for a chain of 
Sc@C82 fullerenes. Empty circles refer to calculated results 
and are fitted using exponential decay law (lines) of the 
form J e- MR - Ro \ where J = 4.0 meV, A = 4.16 A -1 , 
Ro = 3.35 A. Dashed (red online) line refers to J obtained 
from the DFT energy difference between FM and AF states. 
Solid (blue online) line refers to J e s obtained from the Heisen- 
berg model 4t 2 /U e s, where U e s is fitted to be 0.412 eV from 
the value of J at Ro. Arrows show discrete values of R at 
which J is calculated for peapods. Inset shows electron den- 
sity contributed by highest occupied states of Sc@C82 chain. 
Range of isovalues is 0-0.002 e/A 3 . 



energy difference between FM and AF configurations, is 
3 meV per cell (containing two spins). The behavior of 
J as a function of the inter-fullerene separation, R, in 
a Sc@Cs2 chain is plotted in FigJ3] Values of J calcu- 
lated for the peapods at discrete values of R, as indi- 
cated in Fig. [3l coincide with those obtained for the 
Sc@Cg2 chain within the accuracy of the present calcula- 
tions. At these separations the inter-molecular spin inter- 
action is therefore via direct exchange between fullerenes, 
with a negligible contribution from interactions via the 
nanotube. This inter-molecular coupling is much larger 
than the classical dipole coupling of N@C 60 [HI and even 
larger than that computed for defective fullerenes with 
inter-cage links 24J, ensuring > 10 3 two-qubit gate op- 
erations within the decohcrcnce time. This surprising 
result follows from the HOMOs in the Sc@Cs2 chain be- 



ing very extended as illustrated in Fig. [31 The p z orbitals 
on the closest C atoms belonging to adjacent molecules 
overlap in a cr-type fashion. This implies that the ex- 
change interaction could be tuned by varying the sepa- 
ration between the fullerenes in peapods, for example by 
using functionalized fullerenes [25| . Such high J values 
are consistent with recent magnetic susceptibility experi- 
ments on Sc@Cg2 solids [2f|, which show AF Curie- Weiss 
temperature ~ 300 K, consistent with J = 17 meV at 
R = 3.08 A in Fig. 3, which corresponds to the mea- 




FIG. 4: Spin-polarized band structures of the (a) 
Sc@C 82 @(14,7) and (b) Sc@C 82 @(ll,ll) peapods for the AF 
configuration. Ef is the Fermi energy. 



A generic Hubbard- Anderson model, going beyond 
DFT and capable of describing the low-energy spin inter- 
actions, is conjectured and the energy parameters of this 
model are estimated by direct comparison of its mean- 
field solutions with the DFT results. The Hubbard- 
Anderson model may be used directly to estimate the 
Heisenberg exchange between spins along the fullerene 
chain and the Kondo exchange interaction between lo- 
calized spins on the fullerenes and spins of propagating 
electrons or holes in the nanotube. The Hamiltonian is 
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is a creation operator for an electron in the 
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nanotube with band I with quasi-momentum fc, 
a creation operator for an electron in the highest occu- 
pied fullerene band, 7^ is the hybridization parameter 
between fullerenes and nanotube, U is the intra-fullercnc 
Coulomb repulsion energy and m is the number operator 
for the site % of the Sc@Cs2 chain. 

The AF band structures of the peapods are plotted in 
Fig. 4. They display a very narrow but almost perfect 
cosine form for the highest occupied band [l5|. This is 
well described by a simple tight-binding model with the 
dispersion energy Eh = Eq — 2icosfc. It may be shown 
from an AF mean-field solution of the Hubbard model, 
that this dispersion corresponds to E k = Eq — 2£cosfc in 
Eq. (1), where t = VtU, and U ~ 0.53 eV is obtained 
from the energy gap between centres of the highest oc- 
cupied and lowest unoccupied bands in Fig. 4. This so- 
lution is also equivalent to a Stoner model for exchange 
I = U, where J = At 2 /U The behavior of the 

computed hopping parameter, t, as a function of R, in a 
Sc@C82 chain is fitted by an exponential decay law of the 
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form t e- x< - R - Ro \ where t = 20.3 mcV, A = 2.08 A -1 , 
Rq = 3.35 A, showing t 2 and J scaling in the same way. 
For the values of R in the (14,7) and (11,11) peapods, t 
= 17 meV and 2 meV, respectively. 

An effective U is fitted by fixing J = 4t 2 /U e g within 
the Heisenberg model at Rq, giving U c g = 0.412 eV, in 
good agreement with the DFT calculations. Fig. 3 shows 
that J e ff and J are indistinguishable, implying approx- 
imately constant U over the range of R considered. U 
deduced from the DFT band gap is consistent with the 
mean-field solution of the Hubbard- Anderson model and 
the total energy difference between FM and AF solutions. 
With t <C U for the peapods considered, we thus expect 
strongly correlated electron effects with well-defined spin- 
qubits along the Sc@Cs2 chain of the peapods and weak 
charge fluctuations. Therefore, the system is well char- 
acterized by an antifcrromagnetic Heisenberg model. 

In the metallic peapod, we may also estimate the hy- 
bridization coupling parameters 7zfc in Eq. (1) from the 
anticrossing gaps associated with the nanotube bands 
and a narrow fullcrene band. For example, the anticross- 
ing shown in the inset in Fig. |4] (b) yields a coupling 
energy jik=AE/2=5 meV, where AE is the energy gap. 

The weak interaction between electrons on the 
fullerenes and conduction electrons (or holes) in the nan- 
otube, characterized by 7^ discussed above, will give rise 
to Kondo-likc coupling between spins on the fullerenes 
and spins in the metallic nanotube. The energy scale for 
these couplings is given by Jk ~ where E t is a 

charge transfer energy gap, i.e. E t = Eq + U — Ep or 
Ep — Eq [28j|. This gives a typical Jk ~ 0.1 mcV for the 
(11,11) peapod. 

DFT calculations cannot resolve the difference in J ob- 
tained for a Sc@C82 chain and the corresponding peapod 
structure even at large inter-fullerene separation. This 
is consistent with direct exchange dominating RKKY in- 
teractions. Charge fluctuations could be increased by 
cither enhancing the hybridization interaction through a 
decrease in the fullcrene-nanotubc separation, or by tun- 
ing the Fermi energy in metallic nanotubes to approach 
the mixed valence regime in which the charge transfer en- 
ergy tends to zero. This would enhance both the Kondo 
coupling and the RKKY interaction, in competition with 
the direct Heisenberg exchange. 

In conlusion, the perfect doping mechanism of 
Sc@C82 leads to well-defined spin-qubits on the Cs2 cage 
in the peapods considered, coupled via antifcrromagnetic 
Heisenberg exchange interaction. For the semiconduct- 
ing case, the upper and lower Hubbard bands of the 
fullerene chain are little affected by the nanotube and 
occur entirely within the band gap of the nanotube, al- 
lowing excitations of the Sc@Cs2 chain independently of 
the nanotube. Remarkably, in the use of spin peapods 
for quantum technologies, the main function of the nan- 
otubes will be to give mechanical support for the endohc- 
dral fullerenes and to protect them from the enviromcnt, 



rather than to provide controlled interactions between 
the spins. An endohcdral fullerene peapod thus provides 
a candidate nanostructure for spin-chain quantum com- 
puting 0,0. 
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